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I. Introduction 


Many investigations of electro-fluid mechanical interactions consider only 
surface electrical forces, such as is the case of an interface supporting a 
surface charge distribution or of an interface between different dielectrics 
[1], I2J - For these cases, complications resulting from nonuniform equilibrium 
electric field intensities and interfacial curvature occur only for non-planar 
systems. In contrast, recent work has Investigated electrohydrodynamic coupling 
of perfectly insulating fluids supporting volume space charge distributions 
[3], [4]. Such analysis is appropriate if the fluid is highly insulating and 
the dynamics of interest occur over time scales which are short compared with 
those required for initially injected charges to accumulate appreciably at a 
fluid surface. For this case, the volume space charge imposes a nonuniform 
electric field distribution even in planar geometry, which must be taken into 
account when the interfaces move as the interfacial equilibrium electric fields 
which act on the interfaces also change. 

This work examined the propagation and instability characteristics of 
small signal electro-fluid mechanical space charge and polarization waves 
for electrohydrodynamic configurations similar to those of the classic fluid 
problems of the Rayleigh-Taylor instability for superposed charged planar 
layers, and of interfacial capillary oscillations of charged liquid cylinders 
(or jets) and charged spherical drops [5]« A systematic approach was developed 
to handle multi-in ter facial systems of incompressible, inviscid, and perfectly 
insulating fluids in planar, cylindrical, and spherical geometry through the 
use of a general set of relations for perturbation field and flow variables 
on the perturbed surfaces of fluid layers having constant properties [33 • 
Although the methods developed are valid for any geometry , we limit ourselves 
here to systems initially in spherical equilibrium. 

The analysis showed that the electrohydrodynamic coupling for uniformly 
charged layers could be represented as a purely surface coupled interaction, 
even though volume Coulomb forces are present. The pertinent electrical 
equation necessary in this development is Poisson* s equation relating the 
electrical potential <p to the charge density q and permittivity £ (assumed 
constant) 


- 


ra. 

e 


( 1 ) 


The electrohydrodynamic coupling occurs through the Coulomb force density 
f - -qV<f> (2) 
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the methodical approach developed can be used for any system whose force density 
is curl free within a homogeneous layer (q * constant) with the resulting 
coupling occurring through the interfaces as the force density can then be 
lumped with the hydrodynamic pressure. 

this work, will also use this duality with the directly analogous equations 
to (1) and (2) for self •gravitating systems 

V*V * 4irGp (3) 

f - -pW (4) 


where V is the gravitational potential, p is the mass density and 6 is the 
universal gravitational constant (G * 6.67 x 10*11 nt-m^/kg ^) 9 to deduce the 
dynamics and instability characteristics of spherical self-gravitating 
geometries [6]« It is a simple matter to obtain results for self-gravitating 
systems for those electrohydrodynamic problems already solved by making the 
simple substitutions 


+ -► V ; q * p 



( 5 ) 


However, an important distinction between self-gravitating and electrohydro- 
dynamic systems remains. Whereas the charge density can be either positive or 
negative, the mass density is always positive, thus there is always a force of 
attraction between masses, while because of the minus sign difference in 
Eqs. (1) and (3) , like charges repel and opposite charges attract. Because of 
the absence of "negative mass" there are no gravitational analogs to polarization 
effects due to differences in permittivity or electrical shielding due to the 
presence of electrical conductors. The gravitational constant G is independent 
of material properties. Thus with the differences between electrical and 
gravitational systems in mind,' we can immediately write down the solutions to 
those gravitational problems analogous to already solved electrohydrodynamic 
problems using the conversions of Eq. (5).. 

In this work, we will derive the general electrohydrodynamic "prototype" 
relations for a spherical shall and then derive and contrast the dispersion 
characteristics of perfectly conducting and perfectly insulating charged drops. 
Using duality we will then immediately write down the analogous solutions for 
self-gravitating drops. 
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II. Equations of Motion 

1. General Development 

No matter the geometry, the general equations of motion for an incompressible, 
inviscid, and perfectly insulating charged fluid with mass density pj_ velocity 
v, pressure p, charge density q, permittivity £, and electric field E are 

Conservation of momentum - 


p ^ + Vp « qE - i f • 
Conservation of mass - 


( 6 ) 


H - 0 ; 7 . V - 0 (7) 

Maxwell's equations - 

7 X r - 0 ; ¥ - -V4> , (8) 

J 

Gauss * s Law: V * (£E) * q ; - 0 (9) 

Conservation, of Charge: * o (10) 


In the prototype layer, examples of which are shown in Fig. (X) , the liquid 
is homogeneous such that the equilibrium properties of mass density charge 
density q^, and dielectric constant e&, are constant. Then all quantities are 
assumed to have small perturbations from the equilibrium. Subject to the 
constraint of a homogeneous medium the equilibrium variables must obey the 
time independent form of (6)-(10) and must satisfy the boundary conditions 
between regions. 


Because the fluid layer is homogeneous the charge density, mass density, 
and permittivity remain constant in spite of the fluid motion. This is to be 
expected, because any transport of material into a given region leads to a 
transport of material which has the same properties as that previously occupying 
the given region. This statement only applies to those portions of the fluid 
not swept out by interfacial motions. If a point of interest is adjacent to 
an interface, an excursion of the interface could result in an abrupt change of 
properties. However, if surface deflections are considered at a given instant, 
all properties everywhere between interfaces are uniform. 


Denoting perturbation variables with primes, we take the divergence 
of the linearized form of (6) , to yield the set of perturbation equations 

vV - 0 

( 11 ) 

T^ir' * 0 ; it* » p* + q^<f>' 

Thus regardless of the geometry, the problem reduces to solutions of 
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Laplace’s equation, both for the perturbation potential and for the modified 
pressure tt* . . 

2. Generalized Halations For a Spherical Prototype Layer 

We consider the lower prototype layer shown in Fig* (1) with the under- 
standing that the picture is a cross section of a spherical shell* In 
equilibrium the uniformly charged fluid extends over the range 8 < r <. a. 

We denote all variables at the inner interface with a superscript 8 and all 
variables at the outer Interface with a superscript a. The radial perturbation 
displacement's are £$ and When these displacements are zero the equilibrium 
distributions are 

P + q A m constant 
r o 0 


-q.r constant- constant,, 

* J£_+ l + 2 

6e. r 


_ [V , coa8tanC ll 

L*A r 2 J * 


All perturbation variables are assumed of the form 


( 12 ) 


if* m Rat(r)P®(cos i|/)expj [tut - m6] ; m > 0 » a > 0 * m < a (13) 

where P®(cos i|i) are Legendre functions depending on the azimuthal angle ip. 
Substituting into Eq* (11) yields 

#(r) m A^r 11 + A^r - ^ 11 + ^ (14) 

where A^ and A 2 are constants to be determined from the boundary conditions 

$ r (B) - jo£ 6 

05) 

$ r <a) - JuT 

From Eq« (6) we know, that the perturbation velocity is relating to ft(r) as 

$ r (r) - - 1 dft (16) 

dr 

so that the interfacial displacements are related to the parameters A^ and A^ 
as 
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Y 

1 

m 2 

m*" 1 

- (n + l)cf< a + 2) 


*1 

(17) 


pa* 

p6 n " 1 

-<n + l)fr (n + 2> - 


' A 2' 



For the purposes of our analysis , the inverse relations are needed* Once A. 
and are determined In terms of v 1 and v*, substitution into Eq. (14) 
yields the relationship between the interfacial modified pressures and the 
Interfacial displacements as 


where 


ft (a) 

m 

F(o,B) G(a,B) 

ft(B) 


G(B,a) F(B,a) 



n 2 n 

-rzf ' 

(*) £_ + (Z) _s_ 
V n V n + 1 


* n a 1 

$ *] 


rn 


(« 2 P a ) 


G(x,y) 


■(2n + l)j‘ 


n(n + 1) j<j) x - <*) y 


<“ 2 V 


(18) 


(19) 


Note that the analysis implicitly assumes the interfacial displacements 
to be small as Eqs. (15) and (16) were evaluated at the equilibrium positions' 

(a, 6) rather than at the interfaces themselves (a ■+ 3 + £^) * Fortunately, 

because the velocity itself is a perturbation, the difference between evaluating 
it at the interface or at the equilibrium position is second order in the 
perturbation amplitudes. This illustrates the general approach used in linearized 
surface deformation problems. The boundary condition at the moving interface 
is replaced by one at the equilibrium position of the boundary, thus greatly 
simplifying the analysis* 


The analysis is still not complete for as the interfaces deform, in addition 
to perturbing all variables, the equilibrium quantities acting on the Interfaces 
also change. Thus to compute the total first order change in all variables 
evaluated at the interface* linear changes of equilibrium quantities must be 
included. For example, the total linear changes in the pressures are: 


P 


P(a) 



5 ° 


o) 


S 6 


$(&> 




( 20 ) 
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Similarly 


tor 




- 5(o) +$■ 
5 8 - 5(8) + 


r - $(o) - e° 


(r - o) 


I 8 - 5(8) - 


(r - 8) 


( 21 ) 


Using these definitions we have 


p“ s TT(O) - q A $(o) + E a q A C° 
^ = ft(S) - q A $(0) + E e q A t 6 


dE- 

e A S ? * e A ®r (ct) + e A dT 


r * a 


dB- 


e A S r “ e A S r (S) + e A dT 


|S 


r « e 


( 22 ) 


(23) 


Then using Eqs. (20) - (22) we obtain the generalized mechanical relations 
evaluated at the interfaces 


Y 


F(a,8) 

G(a,B) 


* 



f. 


G(g,cO 

F(8,a) 


where F and G are given in Eq. (19) 



(24) 


Similar operations are performed in the solutions of Laplace 1 s equation 
for the perturbation potential $(r) to yield the electrical relations at the 
interfaces 


Y 

e A*r 


B(a»8) 

C(a,8) 


Y + 

e S 6 
e A r 

* e A 

c(e,o) 

B(B,a) 


t + e¥ 



L 1— 



- 


dE 

0 

v-~ 

dr 

r « a 

^0' 


dr 


_ 

r - B 


( 25 ) 
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where 



(26) 


The general relations of (24) and (25) greatly simplify in various limits. 
As the inner radius goes to zero (B 0) , so that the shell becomes a drop , 
the terminal relations reduce to 


lim 
8 -M) 


e A e r 


2 £a 

.a “ *A a C - V 


-e^n ($° + E a £°) + dE Q 
a dr~ 


r « a 


(27) 


As the outer radius becomes very large, (a -*■ *>) , the terminal relations 
become 


lim 
a *► » 


e A S 3 
A r 


^) 2 p a S C 8 - q A $ 6 

(n + 1) 


e A (n + 1) ($ B + E S | 6 ) 

6 


+ s ifo 

dr 


c e 


r - e 


(28) 


Equations (24) and (25) are useful because they relate the interfacial 
variables of pressure, displacement, electrical potential and electrical 
displacement, which appear in interfacial boundary conditions. The boundary 
conditions for all cases include mass and electric potential continuity, the 
second condition being equivalent to the continuity of the tangential component 
of electric field, as well as an interfacial force balance. For perfectly 
conducting interfaces, the perturbation potential must be zero but the surface 
charge imposes a surface force density. Because there are no electric fields 
inside the perfectly conducting fluids, there are no polarization effects. 

In contrast, perfectly insulating fluids with volume charge have nonzero 
interfacial potential and although free charge can make no contribution to an 
interfacial surface force density, there is a surface polarization force if 
the permittivity of the fluid and its surroundings differ. For non-perfect 
conductors (including perfect insulators) no surface charge can be allowed on 



any Interface for an inviscid analysis to be meaningful* as electrical shear 
stresses would accompany such a surface charge density, and our model would 
then necessarily have to Include viscosity (or some other mechanism) to balance 
this shear stress [7], [8]. This requires the normal component of electric 
displacement to be continuous across an Interface for non-perfect conductors* 

An lnvlsdd formulation Is allowed only If the Interface has no surface charge 
or If the Interface Is perfectly conducting, as then the electric field terminates 
perpendicular to the Interface, resulting in no electrical shear force. However, 
the same systematic techniques have been applied to a viscous prototype layer 
in planar geometry so that these limitations may be removed [9]. In general, 
considering viscous fluids greatly Increases the mathematical complexity 
of the analysis but yet has no effect on the conditions for Instability [10]. 
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III . Stability and Dynamics of Charged Spherical Drops 

The modes of oscillation and instability of charged, initially spherical 
drops are obtained for small perturbations by using the generalized relations 
of (24) and (25). Rayleigh's limit determines the maximum amount of charge 
distributed around the surface of an Isolated perfectly conducting drop for 
it to be stable [11]. The Isolated uniformly charged, perfectly Insulating 
drop has historical significance through the liquid drop model for fission 
of the nucleus proposed by Bohr and Wheeler [12] . In this section we wish 
to generalize both these classic problems by immersing the charged spherical 
drops within a uniformly charged region. Figure 2 describes the geometry with 
the understanding that the picture is the cross section of a spherical system. 

In spherical geometry, all perturbations due to small signal interfaclal 
motions are assumed of the form 


£ « Re 5 exp[j(a)t - md) }P m (cos $) ; m>0,n>0,m<n (29) 

n 


For spherical drops we only need to know these relations in the two limits 
where the inner radius tends to zero, 8 + 0, and when the outer radius gets 
very large, a «*, which are described by Eqs* (27) and (28). 


1. Perfectly Conducting Drop With Surface Charge Within a Uniformly 
Charged Region 


We first consider a configuration similar to that in Fig. (2) with a 
perfectly conducting spherical drop of radius R, mass density p, and surface 
tension y with uniformly distributed surface charge density cr f immersed within 
a perfectly insulating charged fluid of infinite extent with charge density 
$2 a 11 * 1 mass density p£. Hie equilibrium electric field distribution is 


E 

r 



q 2 (r - R 3 /r 2 ) 


3e„ 


r > R 


(30) 


so that 


E 


1 


0 





r « R 


«2 


(31) 


Because the drop is perfectly conducting there are no electric fields inside 
(E^ * 0). The perturbation interfacial boundary conditions are 


% * t 2 2 l , \ “$ 2 =0 

P x - P 2 + e 2 E 2 S r2 " h V n " 1} (n + 2) ? “ 0 


( 32 ) 
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From Eqs. (27) and (28) we obtain 


P x - PjfAt/n 

P 2 - -p 2 ia 2 Rg/ta~yi:T'~' (33) 


e 2 S r2 * e 2 B 2^ n + 1) ^ R + s 2 **2 

dr" 


| 

R 


Substitution of Eq. (33) into Eq. (32) yields the dispersion equation as 


2 

We see that the self- fie Id term proportional to Of due to the interaction 
of the surface charge with its own field is always destabilizing, while the 
imposed field term proportional to q*Of due to the interaction of the space 
charge q^ with the field due to Of can be stabilizing if q 2 and a ^ are of 
opposite sign. If O^ » 0, there is no electromechanical coupling; 

space charge in the 


(35) 


We see that the n * 0 mode is not allowed if Pi f* 0 due to the incompressibility 
of the drop, and that the n * 1 mode is neutrally stable. As O^ is increased, ^ 

the first mode to become unstable is n * 2 with critical total charge Q T (Q T - 4irira f ) 

Q t - 8irR<e 2 YR) I/2 (36) 


We examine (34) in Rayleigh * s limit where there is no 
outer region (q 2 * 0), so that Eq. (34) can be written as 





(n - l)/n + 2 - ap\ 



- (n - 1) (n + 2) 


(n + 



(34) 


However if the drop has negligible density so that * 0, which is the 
case of a bubble within an infinite fluid with surface charge on the bubble 
fluid interface then the n » o solution is allowed. In the absence of electrical 
forces, surface tension makes this case unstable by acting to collapse the 
bubble. However, with surface charge present, Coulombic repulsion could stabilize 
the system if 

cy f R ^ e 2 Y > 2 (37) 


For finite q^, the n * 1 mode will be stable only if q 2 <?£ < 0- 
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2. Perfectly Insulating Drop with Volume Charge Within a Uniformly 
Charged Region 

We now consider a perfectly insulating spherical drop supporting a constant 
volume charge distribution within a uniformly charged region, as in Fig, (2) . 

We include the possibility of having a point charge of value Q at r - o, and 
give the drop a dielectric constant different from the surrounding region to 
include polarization forces. 

The equilibrium electric field distribution is 


Q q*,* 

2 + — 

4'fTE^r 3£ x 


r < R 


(<u - q 0 ) R qo r 


\ + * * 
' 2 2 

— + — r > R 

4Tre 2 r ^ e 2 r 

*2 

for the interfacial fields and field gradients 

Q q,R dE. 

-Q 

11 2 h > 

4-!TE 1 R 3e 1 dr 

3 t " 

2716^ 3e 1 


r ** R 

Q q^R dEj 

-Q 2q. q. 

+ , 

4re 2 R 3e 2 dr 

3 + + 

27re 2 R J 3e 2 e 2 

’ 

r » R 


The interfacial boundary conditions result in the relations 

A A A A A 

= C . <I>1 - <t> 2 = ^> 


£ l S rl * e 2 S r2 

h ~ h * E 2 E 2 S r2 " £ l E l S rl “ ^ (n -«(» + » 5 - 0 

From the general relations of Eqs. (27) and (28) we obtain 
“ PjU^RS/n - q 1 $ 


(38) 


(39) 


(40) 


P 2 - -P 2 0) 2 RC/(n + 1) - q 2 $ 


( 41 ) 
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e x e rl - -(e x a/l)($ + E^) + e x dE 1 

dr* 


r « & 


s 2 S r2 * € 2 (n * 15 + E 2 5) + H dE 2 

R dr* 


(42) 


Substitution of Eqs. (41) and (42) into Eq. (40) yields the formidable 
dispersion relation 


P»¥l , p 2 ) 

\ Y /V (n + 1 )/ 


{n - 1) (a + 2) 


< q l - q 2 )R 


Cq l “ q 2 * (2n + 1 ) £ 2 E 2 /R1 


(43) 


y[e 2 + (e x + e 2 )n] 


+ (e 2 - e 

y[e 2 + (e 


X )R 2 e 2 [q^a + b 

e x + s 2 )n] e L 


+ <?2 a " E i^ £ 2 + ^ e l + e 2^ n ^ * a ^ a * 
F 


- v| 


In.Eq. (43), the last term on the right represents the polarization surface 
force due to the difference in permittivities and the fields from the space 
and point charges, the other electric term represents Coulomb forces acting 
even in the absence of a polarization. For simplicity, we examine Eq. (43) in 
various limits. 

BOHR LIMIT (e x * e 2 = £, q 2 ■ 0, Q ■ 0), 


In the liquid drop model of fission of the nucleus, Bohr and Wheeler 
considered an isolated uniformly charged spherical drop under surface tension 
P-2]. Sufficient charging of the drop results in instability, which in the 
context of the nucleus amounts to nuclear fission. Under these conditions 
Eq. (43) can be reduced to 


(44) 


The n * o mode is not allowed for finite due to the drop incompressibility. 
The n * 1 mode is neutrally stable, while the first mode to become unstable as 
q^ is increased occurs for n * 2 with total charge Q^(Q T * 4irR3q 3 /3) 



/»! + »2 ' 

1 - (a - 1) 

(2n + l)(a + 2) - 2q 2 R 3 


\a (a + 1)/ 

1 (2n + 1) 

3ey . 


Q x - 8 ttR(5£yR/6) 


1/2 


(45) 


We see from Eq. (45) compared to Eq. (36) that a perfectly conducting drop 
could support slightly more total charge before becoming unstable than the 
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perfectly insulating drop* 

NO POLARIZATION EFFECTS (^ *~e 2 5 £) 

There are no polarization forces if the dielectric constants of the two 
regions are title same. Generalizing the Bohr limit by allowing finite q 2 and 
Q simplifies Eq. (43) to 

A P 2 \ /Y \(n - IHn + 2) + (q x - q 2 ) 2 

“\n + (n + 1)/ * ( R 3j e(2a + 1) 



£ 


The second term on the right of Eq. (46) is a self -field term and is 
always stabilizing while the last term is an imposed field term and can be 
either stabilizing or destabilizing depending on the relative charge polarities. 
A sufficient condition for stability is for the last term in (46) to be negative. 

NEUTRAL DROP (qj4irR 3 /3 + Q * 0) 

If the total charge within the spherical drop is zero, the equilibrium 
inter facial electric field is also zero. Under this condition even if the 
permittivities of the two regions are different* there would be no polarization 
effects as the interfacial electric field is zero. However, due to the electric 
field gradient, there is still an electromechanical coupling so that Eq. (46) 
is still appropriate with the last term being zero. This system is stabilized 
by the space charge. In fact, for the case of a bubble where Pi m 0 the 
destabilizing nature of surface tension for the n «* o mode can be opposed so 
that the system becomes stable if the difference of charge densities are of 
sufficient magnitude such that 

<q x - q 2 > 2 > 2SY/R 3 (47) 

UNIFORM SPACE CHARGE DENSITY THROUGHOUT (q 1 - q 2 ) 

The interplay of the polarizability and the space charge is demonstrated 
by considering a case in which the respective space charges are equal (q^ * q^) . 
Then the second term on the right of Eq. (43) drops out. If the relative 
permittivities were also equal, there would be no electromechanical effect. 
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IV. Gravitational Stability of Spherical Drops 

For the space charge analysis, the mass density only contributed through 
Inertial terms. For the self-gravitating counterpart, the mass density will 
also contribute via the force law of Eq. (4). We now consider the gravitational 
analog to Fig. 2 where a fluid sphere with mass density p- is immersed within 
an infinite medium of mass density p 2 - A point mass M is at the origin r » 0. 
The equilibrium gravitational field in the negative radial direction at the 
interface is 


8 


MG 


41TGP 1 R 


(48) 


We can use the results obtained for a perfectly insulating drop with 
volume charge within a uniformly charged region by using the dual relations 
of Eq. (5) , substituted into Eq. (43) with to obtain the self-gravitating 

dispersion relation 



(49) 


KELVIN’S MODES 

Kelvin first considered this problem in the limit, M * 0, P 2 * 0, and 
Y * 0 [13] . Then Eq. (49) describes the oscillation frequencies^! an inviscid 
liquid globe under its own gravitational field. 

ti) 2 - 8irGp^a(n - 1) (50) 

3(2n + 1) 


HOLLOW DROP 

In the other extreme, a spherical void (M * 0, Pi * 0) within a large 
medium is always unstable for n ■ 0 and n - 1 as the right hand side of Eq. (49) 
is always negative and thus unstable. Unlike the space charge dual, self- 
gravitating forces can never stabilize a spherical bubble. 
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V. Concluding Remarks 

Because the perturbation electrical and gravitational force was curl free 
within a homogeneous layer, the broad class of interactions involving liquids 
with discrete stratifications in mass density-, charge density, and permittivity 
are representable as surface coupled interactions. The transfer relations derived 
are useful since they relate interfacial variables that appear in the boundary 
conditions. The techniques developed here may also be used for other perturbation 
volume forces which are curl free within a homogeneous layer, or for surface 
forces which only act at an interface. Any system which is described by Poisson’s 
equation and a Coulomb force law can be immediately solved using similar dual 
relations as in Eq. (5). 

Since the fluid's are modeled as inviscid, force equilibrium at the 
interface required only a normal stress balance. No electrical shear stresses 
are allowed because of the absence of viscous shear forces to oppose the 
electrical forces. A more general analysis which would allow electrical shear 
forces must include viscosity. Then the mechanical terminal relations must 
include relations between the shear and normal forces and the shear and normal 
interfacial displacements. Force equilibrium at the interface requires both a 
normal and shear stress balance. This results in the mechanical transfer 
relations becoming 4x4 matrices, rather than the simpler 2x2 matrices 
considered here, greatly increasing the algebraic complexity. 
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Fig. 1 Prototype layers of incompressible, inviscid, perfectly 

insulating fluids, supporting uniformly distributed charge 
densities, (a) planar geometry, (b) cylindrical or spherical 
geometry. 



Fig. 2. Cross section of an Initially spherical drop of radius R 
with charge density mass density P. and permittivity 
placed within an Inf inite medium with charge density <u> 
mass density and permittivity e~. A point charge (Q -X) 
is at r » 0. ■ * 



